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Abstract. The Abelian Sandpile Model is a cellular automaton whose 
discrete dynamics reaches an out-of-equilibrium steady state resem- 
bling avalanches in piles of sand. The fundamental moves defining the 
dynamics are encoded by the toppling rules. The transition monoid 
corresponding to this dynamics in the set of stable configurations is 
abelian, a property which seems at the basis of our understanding of 
the model. By including also antitoppling rules, we introduce and in- 
vestigate a larger monoid, which is not abelian anymore. We prove a 
number of algebraic properties of this monoid, and describe their prac- 
tical implications on the emerging structures of the model. 



1 Introduction 

In 1987 Bak, Tang and Wiesenfeld [1] proposed a simple cellular automaton model of 
sandpile growth (since now BTW model) as an example of self-organized criticality 
(SOC). These systems should be characterized by a dynamics which sometimes, in 
an apparently unpredictable way, shows bursts of activity, avalanches in the sandpile, 
which eventually drive the system into an out-of-equilibrium steady state. The state 
is a steady state, in the sense that overall properties stay unchanged during time 
evolution, and out-of-equilibrium because these systems are open and dissipative and, 
therefore, they require input from outside at a constant rate to balance the dissipation, 
and present a spatial distribution of current. These steady states are critical as they 
exhibit long-range correlations with power-law decay. 

Consider a rectangular portion of a square lattice of size L^xLy, and let n = L^Ly. 
To each site i — (x, y) is associated a non-negative integer height variable (sometimes 
called mass) Zi. A configuration is the collection of all heights and will be denoted as 
a bold letter: z = {zi}"^!. If all Zj's are smaller than 4, the configuration is stable. 
Otherwise, a procedure called relaxation gives univocally a new configuration TZ(z). 
The relaxation is performed through a sequence of elementary moves, called topplings: 
if Z(^x y-) > 4 (in this case we say that the site {x,y) is unstable), then we decrease 
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^{x,y) by 4, and increase by 1 all the values at neighbouring sites, Z(^x±i.y) f^nd z(^x y±i^ 
(nothing is done at those coordinates (a; ± l,y), {x,y ± 1) falling outside the grid). 
Topplings are performed at the various unstable sites, until when the configuration 
obtained is stable. 

It is remarkable, and important, that the final result does not depend on the order 
in which the topplings are performed. In particular, if i and j are unstable sites in the 
configuration z, then j is still unstable after a toppling has been performed at i (and 
viceversa). Indeed, all the valid sequences of topplings differ only for their ordering 
(although not all the permutations of a valid sequence are valid sequences). 

Note that, if we call |z| = X]i=i ■^i ^^e total mass of the configuration, the mass 
decreases for topplings occurring at the boundary of the grid, and is conserved for 
topplings occurring in the bulk. 

Call Qi the operator which maps a stable configuration to a new one by the addi- 
tion of a grain of sand at the site i, followed, if it becomes unstable, by a relaxation. 
We have n operators a^'s, acting as a transition monoid on the set of stable configu- 
rations. As first pointed out by Dhar [2], these operators commute, i.e., for all z, the 
configurations aiQjZ and OjOiZ coincide. As a result, the transition monoid is abelian. 
Under the Markov Chain obtained by acting with the a^'s, chosen at random, the 
sandpile is ergodic in a subset of all possible stable configurations. According to the 
general notion used in the analysis of Markov Chains, these configurations are called 
recurrent, and they are characterized by the absence of forbidden sub- configurations 
(FSC). All recurrent configurations occur with uniform probabilities in the steady 
state of this dynamics. A stable configuration which is not recurrent, i.e. that has 
zero steady-state probability, is called transient. As we will recall later on, the no- 
tions of recurrent and transient configurations can also be posed in purely algebraic 
terms, with no need to refer to probabilistic notions. 

In the steady state the fluxes of sand added to and flowing out of the system must 
balance. This balancing occurs however in a non-trivial way, as the distribution of 
the size of the avalanches, i.e., the number of topplings required in the relaxation 
process at a given step of the Markov Chain, has an algebraic tail, although the 
precise exponent is still unknown [34 5 6. 

The algebraic behaviour of the distribution of avalanche sizes suggests a relation 
with a state at criticality. Indeed, there is an equivalence between the sandpile model 
and the q ^ limit of the g-state Potts model, at criticality, through a common com- 
binatorial description in terms of (rooted) spanning trees [7] . Therefore this sandpile 
critical behaviour is in the same universality class of a conformal field theory (CFT) 
with central charge c = —2, a relation which has allowed to obtain many exact re- 
sults on the critical exponents associated to quantities that, differently from avalanche 
sizes, are more easily related to fields in the CFT. For example, the correlation of 
height variables at two sites decays only algebraically with the Euclidean distance 
between them [5], the distribution of the heights has corrections, depending on the 
Euclidean distance from the boundary, that decay only algebraically, and, consistently 
with the fact that the CFT at c = —2 is logarithmic, certain two-point correlation 
functions, or one-point boundary correlation functions, show the logarithmic correc- 
tions to power-law decays that one expects from the presence of Jordan blocks in the 
representation of the operator Lq of the Virasoro Algebra (see |9I10I11I12I13I14|T5] 
for a collection of both rigorous and non-rigorous theoretical results in this direction) . 

The dynamics of the sandpile, historically introduced first in [1] for the square 
lattice, is so simple and natural that it has been soon generalized to an arbitrary 
graph P], an extension that allows to elucidate the main algebraic elements of the 
theory, and also independently studied by combinatorialists, under the name of Chip 
Firing Game |16I17I18I19I20I21| , in relation to the Tutte polynomial of the graph (not 
surprisingly, as the Tutte polynomial of a graph is a reformulation of the Potts Model 
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partition function, the adjacency between variables being encoded by the graph, see 
e.g. [^). In particular, a series of papers |19I20I21] shows how the generating function 
of recurrent configurations, weighted according to their total mass, gives access to 
larger class of enumerating problems, in particular, of spanning connected subgraphs 
and spanning forests. 

Also our personal interest was induced by the relation of the sandpile model on 
an arbitrary graph and the g — >■ Potts model, because we have introduced a novel 
representation in terms of fermionic variables of a model of random forests |23l24j . 
that generalizes the Kirchhoff Matrix- Tree Theorem, corresponding to the case of 
spanning trees. 

More recently, we have been strongly impressed by some general results which 
have been obtained for dynamics where all the sand is deterministically injected in 
a limited region j25l26j . producing complex and beautiful patterns which display 
allometry. This has led to a long-term research project, of which some first results are 
reported in [27) . We believe that an important ingredient to better understand the 
emergence of these patterns, even for the ordinary Abelian Sandpile described above, 
is to understand two distinct, and apparently unrelated, generalisations of the model. 
The first more general toppling rules concern multitopplings, i.e. topplings associated 
to clusters of sites [28129] . The second one concerns antitopplings, where instead of 
adding sand, and then relaxing, the sand is removed, and then an inverse relaxation is 
performed. An analogous of this latter elementary move had already been considered 
in [30] ■ with the aim of determine theoretically certain critical exponents associated 
to the model in two dimensions, and, in the same variant considered here, has been 
studied later on, mostly numerically, in [31] . 

In this paper we will mostly present theoretical results for the transition monoid 
associated to the sandpile in which both topplings and antitopplings are allowed 
(remark that this monoid is no more abelian). These results are in the form of non- 
trivial algebraic relations (see e.g. Theorem [1] later on). Along the paper we will see 
how one of the theorems in the theory of antitopplings (Theorem [5]) has a relatively 
simple proof only in the framework of the theory of multitopplings, this explaining 
the need of treating the two generalisations simultaneously. 

The paper is organised as follows. Section [2j by far the longest of the paper, gives 
in Section \TJ\ a first reminder of the properties of the Abelian Sandpile Model, then 
summarizes in Section a collection of new theoretical results, finally in Section 1^31 
more subtle aspects of the Abelian Sandpile Model are reviewed in detail. Due to a 
natural involution symmetry between toppling and antitopplings, the latter are intro- 
duced from the very beginning, even in the review parts. Conversely, the introduction 
and study of multitopplings is postponed to Section |3l where Theorem [2l presented 
in Section 12.21 is proven. Section U presents numerical examples (on the BTW set- 
ting) that motivate the study of several sandpile dynamics involving simultaneously 
topplings and antitopplings. 



2 Algebraic formalism 

The reader interested in the general theory of Abelian Sandpiles will find in [52] 
a beautiful review. An updated version is in [33] . and other details are reported 
in [33]. In Section [2?T] we shall first report the basic facts that we need for our new 
developments. Note that, differently from the treatments in the references above, we 
introduce antitopplings, as well as several other "anti-" quantities, from the very 
beginning, and we choose a notation that highlights the role of a simple natural 
involution, that exchanges ordinary and the conjugated "anti-" quantities. 
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Facts already described in [32] are reported here for both classes of quantities. 
Indeed, because of the involution symmetry, the action of the conjugated operators 
alone cannot lead to any special novelty. Nonetheless, non-trivial aspects will arise 
when ordinary and conjugated operators simultaneously enter the game. The first 
simple properties in this wider context are discussed along Section 12.11 while the 
most remarkable facts are anticipated in Section 12.21 and proven in the body of the 
paper. Section [2731 presents further aspects of the theory as developed in [32], that 
are not required for the understating of the statements of our results, but are useful 
in their proofs. 

2.1 Definitions and basic facts 

Let n be an integer, the size of the system. It is often useful to think at the system as 
a graph with n sites, and the set of toppling rules in terms of the adjacency structure 
of the graph, thus we will call sites the indices i G [n] = {1, . . . , n}. Consider vectors 
z £ Z", where we will use the partial ordering ^ such that u ^ v if < for each 
i G [n]. We will define the positive cone fl as the subset of w e Z" such that w ^ 0, 
where has vanishing entries for all i. 

An abelian sandpile A = ^(Z\,z, z) is identified by a triple (Z\,z, z), that we now 
describe. The vectors z and z are the collection of upper- and lower-thresholds^ {zi} 
and {Zj} respectively, and are constrained to the condition — > for all i. Define 
the spaces 

^+ = {z G z" I z ^ z} - {8)k„^. + 1, • • •} ; (1) 

i=i 

5_ = {z G Z" I z ^ z} = . . , - 1, zj ; (2) 

i = l 

5 = {z G Z" I z ^ z ^ z} = ^{zi, . . . , = 5"+ n S*- . (3) 

Thus 5+ and —S- are translations of i?, while 5' is a multidimensional interval, and 
has finite cardinality. 

We also have a n x n toppling matrix A, with integer entries, that should satisfy 
'Zi — z^ -\- 1 > An > 0, and Z\y < for i ^ j. We say that the sandpile is tight if 
"Zj — z j + 1 = Aii for all i . Define the spaces 

S+ = (S){zr- Au + l,z^- Au + 2,...}; (4) 

i=l 

S'_ = (g){. . . + Au ~ 2,z,, + Au ~ 1} . (5) 

i=l 

We have 5+3 5^, and S- 3 S'_ in general, while 5+ = S*^ , and S- = S'_, if and 
only if the sandpile is tight. 

We further require dissipativity, that is b~ — Aij > 0. As seen below, b~ is the 
amount of mass that leaves the system after a toppling at i. The requirement that 
the toppling matrix is irreduciblq^ ensures that the avalanches are finite (and that 
detZ\ > 0). For future utility, we also define 6^ = X^i^ii' which is the difference 
between the amount of mass that can leave the site j in a toppling and the mass that 
can be added if all other sites would make a toppling. A site j where 6^ < is said to 

^ This means that for every jo there exists a sequence {jo,ji, ■ ■ ■ ,je) such that Aj^ j^^^ < 
for all < a < ^, and b~ > 0. 



Multiple and inverse topplings in the Abelian Sandpile Model 



5 



be greedy or selfish. Clearly = hf , so also the 6^'s are, on average, higher 

than zero, however positivity on the 6^'s is not implied directly by the positivity of 
the 's. We require the absence of greedy sites as an extra condition, whose utility 
will be clear only in the following. A sandpile is said to be unoriented if Z\ = 
(and thus b~ = b+). 

The above conditions complete the list of constraints characterizing valid triples 
(Z\, z, z). The special case of the BTW sandpile corresponds to A being the discretised 
Laplacian on the square lattice. An = 4, Z\y = —1 if d{i,j) — 1 and Aij = if 
d{i,j) > 1, where d{i,j) is the Euclidean distance between the sites i and j, and 
"Zi = 3, Zj = for all i. 

The matrix A is the collection of the toppling rules, and is conveniently seen as 
a set of row vectors = {Z\ij}i<j<„. Denote by ti the action of a toppling at i. If 
such a toppling occurs, the configuration z is transformed according to 

t,z = z - A, . (6) 

A site i is positively-unstable (or just unstable) if Zi > 'zi. In this case, and only in 
this case, a toppling can be performed at i. Note that, after the toppling, it is still 
^ (more precisely, Zi > li — An), while Zj's for j ^ i have not decreased, thus the 
topplings leave stable both spaces 5*+ and S'^ . The relaxation operator TZ is the map 
from 5+ to S, coinciding with the identity on S, that associates to a configuration 
z G S+ the unique configuration TZ{z) G S resulting from the application of topplings 
at unstable sites. Unicity relies on the abelianity of the toppling rules, i.e. that, if i 
and j are both unstable for z, j is unstable for tiZ, and thus crucially relies on the 
fact that Aij < for i ^ j. 

A site i is negatively-unstable if Zi < z^. In this case, an antitoppling can be 
performed at i, with the rule 

4z = z + A,. (7) 

We use deliberately the symbol instead of because, although the effect of the 
linear transformation ([7]) is just the inverse of the effect of ^ , the if conditions for 
applicability of the two operators are different. 

Now antitopplings leave stable both spaces S- and S'_ . The antirelaxation operator 
is the map from S- to S, coinciding with the identity on S, that associates 
to a configuration z e S'_ the unique configuration 7?.^(z) e S resulting from the 
application of antitopplings at negatively-unstable sites. As a matter of fact, the 
involution 

i : z— s>z + z — z (8) 

exchanges the role of operators with and without the f suffix, i.e., for the operators 
above and all the others introduced later on, we have A^{z) = lA{lz). 

Note that, for a configuration z G Z", we cannot exchange in general the order- 
ing of topplings and antitopplings (i.e., if i and j are respectively positively- and 
negatively-unstable for z, j might be stable for t^z, and i might be stable for i]z). 
Consistently, we do not define any relaxation-like operator from Z" to S, as it would 
be intrinsically ambiguous, and even not finite, even when the toppling matrix is 
dissipative 0. 

Remark, however, that the definition of TZ can be trivially extended in order to 
map unambiguously Z" to S-, by letting it produce a toppling only on unstable 

^ For example, in a sandpile that has all 1 = 2, all z — 0, all diagonal Aaa ~ 3, and 
non-zero Aab with a = i,j,k given by the list {Zi^ , zljfe, Zifc,, Zi,,i ^ , Zijji ^ , Zifcfej ^ } and aU 
equal to —1, in the configuration z depicted below one can repeat infinitely many times the 
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sites, and, analogously, Tt^ to Z" to 5*+ (and both TZ'^TZ and T?.??,^ map Z" to S, but 
7e1'(7e(z)) ^ 7^(7^1■(z)) in general). 

On the space Z" we can of course take linear combinations, and define the sum 
of two configurations, z + w, and the multiplication by scalars fc G Z, kz. For generic 
z and z these operations do not leave stable any of the subsets S± and S. But the 
sum and difference, z + w and z — w, are binary operations in the following spaces: 

z + w : S+xn^S+; S'_^xn^S'+; (9) 

z - w : S-X f2^ S'_ x [2 ^ SL ■ (10) 

Thus, in particular, the two maps ^Z(z, + ^/v) and TV (z. — 'w) can act from S x f] into S. 
This suggests to give a special name and symbol to the simplest family of these binary 
operations, seen as operators on S. Call the canonical basis of Z", i.e. (e,i)j — Sij. 
Define the operators di such that a^z = z + e^, and introduce the operators of sand 
addition and removal 

Oi = TZdi ; a\ = TZ"^ . (11) 

The OiS commute among themselves, i.e., for every z, aiOjii — OjOiZ = 7^(z + ei + ej). 
Similarly, the aj's commute among themselves. More generally, for any z e Z" and 
we/?, 7^(7^(z) + w) = 7^(z + w), and 7^■|■(7^1■(z) - w) = 7^t(z - w), this implying 

n n 

7^(z + w) = (n(a.r)z; 7^t(z-w)= (n(al)"")z- (12) 

1=1 i=l 

This can be seen by induction in z, as 7?.(z + (w' + e.;)) = TZ{TZ{'z, + e^) + w') = 
TZ{{aiz) +w')). For later convenience, for w G i7, we introduce the shortcuts 

n n 

aw = - 7^( • + w) ; al^ Y[{alr = 7^t( • - w) ; (13) 

1=1 1=1 

(note that the order in the products does not matter). 

These properties have an important consequence on the structure of the Markov 
Chain dynamics introduced in Section [TJ at each integer time t a site i{t) is chosen 
at random, and z{t + 1) = ai(t)z(i). As long as we are interested in the configuration 
z(ifin) for a unique final time tfin, for a given initial state z(0), it is not necessary to 
follow the entire evolution z{t), ioi < t < tun, but it is enough to take the vector 
w = ei((), and evaluate TZ{z{0) + w). Thus, the final result of the time evolution 
depends on the set {i(i)}o<t<t,„ax of moves at all times, in a way which is invariant 
under their permutations. 



cycle it]t^tlt/,t^)z = z. 
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Note however that, similarly to topplings with antitopplings, also each operator 
of sand addition and removal do not commute, i.e. a^OjZ ^ flja^z in general, even for 

z £ 5'. Therefore, in a Markov process involving both a^'s and aj's, in order to know 
the final configuration, it is necessary to follow the full trace of the time evolution. 
We will describe and briefly investigate a collection of dynamics of this kind, for the 
sandpile on a square lattice, in Section |4l Furthermore, see |31) for a first extensive 
investigation of a dynamics in this family. 

Beside the commutativity relations a^aj — aja.i (and conjugated ones), there is a 
collection of n relations, encoded by the toppling matrix: for any i, when acting on 
configurations such that Zi >'Zi — An, (respectively, such that Zi < z_^ + An), we have 

«^-n«7^''; (ai)^^^=n(«i)-^-- (14) 

This because tiaf" — Yij^i '^j such configurations, and the site i is certainly 

positively-unstable after the application of af" . As a corollary, relations valid for 
generic configurations in 5+ , and in S- , respectively, are 

o-i =0-1 11 (15) 



which are concisely written as 



2i-2. + l _ 

°i ^ '^-A,+(¥i-z, + l)e. 



(17) 

(«Ir'""' + '=«'-A, + (..-^. + l)e.- (18) 

Remark that the two sets of equations and do coincide if the sandpile is 
tight, i.e. if — z j + 1 = An for all i. 

In the previous paragraphs we only considered sums z + w, in which one of the two 
configurations is taken from a space among 5, S'±, or S'j^, and the other one from i7, 
or other spaces with no dependence from z and z. 

In such a situation, we have a clear covariance of the notations under an overall 
translation of the coordinates. I.e., for any r G Z", under the map z — >■ z + r, we 
have an isomorphism between the sandpile ^(Z\,z, z) and „4(Z\, z + r, z + r). We call 
a gauge invariance of the model the covariance explicitated above, and a gauge fixing 
any special choice of offset vector r. 

As always, a gauge fixing reduces the apparent number of parameters in the model. 
Some special choices simplify the notations in certain contexts. For example, we can 
set z = 0, so that 5*+ = i7, that we call the z = gauge, or Zi = An — 1, so that 
S'^ = f2, that we call the positive-cone gauge. 

In particular, the positive-cone gauge is the most natural one in the Abelian 
Sandpile in which antitoppling are not considered, as in this case the parameters z 
do not play any role. Conversely, the covariant formalism is the only formulation that 
does not break explicitly the involution symmetry implied by ([8]). 



2.2 Statement of results 



Two new theoretical results are extensively required for the analysis of interesting 
dynamics, in Section |4l The first one is the following 
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Theorem 1 For every i G [n], acting on S+, 




(19) 



and, acting on S^, 




(20) 



This theorem will be proven in the next subsection. As an aside, we know that 
Theorem [1] is the tip of the iceberg of a much more general family of identities, with 
operators a,- and a\ replaced by and a]^^ for w G that we plan to elucidate in 
forthcoming works. Note for example how, acting on S+, the identity a^aX^a^ = 
holds for simple reasons if w is recurrent, due to the fact that a^al^a,^z ^ a^'z, 
both sides of the equation are stable recurrent, and there exists a unique stable recur- 
rent representative in each equivalence class (see later Section 12.31 for the pertinent 
definitions). 

An immediate corollary of Theorem [T] is the following 
Corollary 1 For every i G [n], acting on S, a\a^ and a^a\ are idempotents, i.e., 



Indeed, the set S = S+ H S- is left stable by the action of the monoid. It is enough 
to multiply equation (jl9p by a|, from the left or from the right respectively for the 
two claims, or, alternatively, multiply (j20p by , from the right or from the left 
respectively. 

The simplicity of Corollary [T] may suggest that abelianity is restored at the level 
of these idempotent combinations. This is not the case. No pairs of distinct opera- 
tors in the set {ajaj^, . . . , aJjO^, a^^aj, . . . , a„ajj} commute with each other, in general. 
Nonetheless, a few interesting facts are found. 

For a finite set / C [n], call Afj = {aja^j.^i- For X C Z", call Afi[X] the set of 
y G Z" such there exists a configuration x £ X and a finite sequence (ii, . . . ,ik) of 
elements in I such that a] a, a] a,- • • • a] a, x = y, that is Afi\X] is the set of 

''k •'fc •'fc — 1 — 1 ^1 '•1 " ' ^ -' 

possible images of X under the action of products of operators in J\fj. 
The second new theoretical result of this paper is 

Theorem 2 Consider a sandpile A{A, z,z) such that A = A'^ and Aij G {0,-1} 
for i ^ j. For any I C [n], and any z G 5*+, there exists a unique state y(z, J) in 
J\fi[{z}] such that ala^y — y for all i G /. For any state x G A//[{z}], we also have 



Thus, this theorem shows that certain collections of idempotents have a well-charac- 
terised set of common fixed points. The portion of this set accessible from any con- 
figuration z has cardinality exactly 1. 

The statement of this theorem can be translated in terms of a stochastic evolution. 
Consider a Markov Chain in which the initial state is z(0) — z, and at each time t 
an element a\^a,^^ G A// is chosen (with non-zero probabilities for all elements), and 

z{t -I- 1) = a\^a^^z{t). Then the theorem states that this Markov Chain is absorbent, 
on an unique state y, and in particular, no matter the evolution up to some time t, y 
is still accessible from z{t) (and in fact it will be reached at some time). A dynamics 
of this kind will be described in Section 14.31 

The theorem above will be proven in Section [31 Furthermore, the state y(z, /) 
will be shown to have a further characterization, in terms of a multitoppling Abelian 
Sandpile associated to the original system. 




a\a^, and (0^4)2 




y eM[{x}]. 
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2.3 Further aspects of the theory 

There exists a natural equivalence relation on vectors z G Z", that partitions this set 
into det A classes, which are affine subspaces of Z", all isomorphic under translation. 
This notion was first introduced and studied in [35]. We recall here only briefly the 
easiest facts. 

We say that z w if there exists T G Z" such that z — w = TA. In particular, 
as z — tiz — — (z — t|z) = = BiA, we have z ^ ^ 7^^(z). Analogously, 

as af"z - {Hj^iaJ^'') z = -eiA, we get af"z ^ {H^^iaJ'^") z for all i, as 
it should at the light of (fT4)) . The dissipativity condition on the toppling matrix 
ensure that detZ\ ^ 0. Thus A~^ exists, and it is evident that z ~ w if and only if 
zZ\-i - wZ\-i e Z". So, the fractional parts Qf"'''\z) = (zZ\-i), - [(zZ\-i),J , called 
the charges of the configuration z, completely identify the equivalence class of z. As 
a corollary, Qf'^'^'^\z) = Qf"^'^\t^z) = Qf'^^'^\tlz) (when t^ or 4 are applicable to z). 
For future convenience, we also define Q(z) = {zA~^), so that Q.f'^'*'^^ = Qi — [Qi\ 0- 

The set S of stable configurations is divided into the two subsets of stable transient 
and stable recurrent configurations, S = T U R. Several equivalent characterizations 
of recurrency exist, some of which extend naturally to 5+, and even to the full space 
Z" (this is also our choice). 

In particular, we give three definitions, all valid in Z". Under all definitions, a 
configuration is transient if it is not recurrent. 

Definition 1 A configuration z is recurrent by identity test if there exists a permu- 
tation a e 6n such that ^^(n) ■ • ■ ^(t(2)^(t(i)(z + b+) ~ z is a valid toppling sequence. 

Definition 2 A configuration z is recurrent by toppling covering if there exists a 
configuration u, such that z = ti^^ ■ ■ - ti-^u is a valid toppling sequence, and at least 
one toppling is performed at each site. 

Definition 3 A configuration z is recurrent by absence of FSC's if for every set I 
of sites, there exists i £ I such that Zi >'Zi — X^je/ ^i*- 

The reasonings of the following paragraphs will prove, among other things, that these 
three definitions are equivalent. 

First of all, a configuration recurrent by identity test is also recurrent by toppling 
covering (one can take u = z + b+). 

Note that, if u = 7?.(v) and v is recurrent by toppling covering, also u is recurrent 
by toppling covering. Furthermore, as after a toppling ti one has Zi > Zi — An, all the 

^ This topic would deserve some more lines. The charges Q^l^^'^^ (z) have two disadvantages: 
they are redundant in general, and require "fractional parts", instead of integer arithmetics. 
Both these issues are solved through the identification of a different set of charges, Ii{z). A 
result of the classical theory of the Smith Normal Form is that any non-singular n x n matrix 
A can be decomposed as Z\ = ADB, where A and B are square matrices with determinant 
±1, and D = diag(di), with di > d2 > ■ ■ ■ d„, all matrices having integer coefficients. The 
decomposition is not unique, but the di's are, and have a further characterization that makes 
unicity manifest. Call gi the greatest common divisor among the minors of A of size n — i, 
and pose g„ = 1. Then di = gi-^/gt. One defines /f'''(z) — {A~^z)i, and < li < di, 
with li = Jj-"' modulo di. Dropping the trivial charges, associated to di — 1, the new set of 
charges is non-redundant and completely identifies the equivalence classes. The new charges 
are related to the previous ones through I — DB Q. A more detailed treatment can be found 
in the original reference, We do not make use of these improved charges in the present 
paper. 
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recurrent configurations, either by toppling covering or by identity test, are in fact 
contained in S'^. This last reasoning is a first example of forbidden sub- configuration 
(FSC), whose generalisation involves more than one site at a time. For z e Z", and 
/ C [n], define z|/ as the restriction to the components Zi with index « G /. 
For a set /, define the vector fmax(-f) G as 



(21) 



We say that the pair (/, z|/) is a forbidden sub- configuration for z if z ^ fmax(-^)- It 
is straightforward to recognize that z is recurrent for absence of FSC's if and only if, 
for all /, zj/ 2^ fmax(J'), thus legitimating the terminology. 

A collection of the pairs (/, fmax(-^)) with smallest |/| in the BTW sandpile is as 
follows: 
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The connection between the definitions of recurrent by toppling covering and by 
absence of FSC's is given by the following statement, slightly more general that what 
would suffice at this purpose. 



and V = ti^ti^_^ 



Proposition 1 Let \i ^ S\ 

for all B such that \B ^ A\ < I, v|b ^ fmax(S). 



•ijjU. Define A = Ui<a<fe{*i}- ^^^'^ 



This proposition implies as a corollary (for \B \ A\ — 0) that configurations which 
are recurrent by toppling covering are also recurrent by absence of FSC's. The case 
\B \ A\ = 1 also emerges naturally from the proof. 

Proof. The claim v\b 2^ fmax(S) can be restated as the existence of s G _B such that 
^'s > fmax{B)s- We will producc a valid choice for s. 

If |i? \ A| = 1, choose as s the only site in B and not in A. Let u' = u and 
r(s) — in this case. Otherwise, for all i € B, call t(«) the maximum 1 < a < fc such 
that ia = i, then choose as s the index realising the minimum of r(i), i.e., the site of 
B that has performed its last toppling more far in the past. Call u' = ii^,,) • • -iijU, 
the configuration obtained after the last toppling in s. 

Note that, as u £ 5^ and this space is stable under topplings, u'^ >'Zs ~ ^ss + 1- 
In the remaining part of the avalanche, no more topplings occur at s. Furthermore, 
all the other sites j € B do topple at least once, at a = T(j). Thus we have 



k 

a=T(s) + l 



^iaS ^ Zg ^SS 



k 

a=T{s) + l 



1 



E 



r(s)<a</c 



(22) 



□ 



The comparison with the definition (pij) of fmax(S) allows to conclude. 

We have now the ingredients to prove Theorem [T] 

Proof of Theorem [T] The two equations are related by the involution, so we only 
prove (fTO|) . I.e., for all z G 5*+, a^ala^z = a^z. First of all, as, for z G 5*+ and w G i7, 
= avi,TZ{z), we can restrict our attention to z G S*. If Zi < we have aja^z — z. 
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and our relation follows. If Zi — 'Zi^ the avalanche due to the action of performs 
at least one toppling at i. By Proposition [H y = a^z has no FSC's with / = {i} 
or / = That is, either y.i > Zi — An, or — Zi — An and, for all j ^ i, 

Uj ^ '^j ~ '^jj ^ ^ij- By direct inspection, in the first case a^ajy — aiCi^^y ~ y and 
in the second case a^a[y = t^dit\a^ y = y. In both cases, our relation follows. □ 

For a sandpile A(A,z,z), call A\i the sandpile described by the toppling matrix Ajj 
the principal minor of A corresponding to the rows and columns in /, and by the 
threshold vectors z|/ and z\j. Also call b+(/) the vector b+ associated to Aj j. The 
observation that ultimately allows to relate the definitions of recurrency by absence 
of FSC's and by identity test is the fact that 

fmax(/)+b+(/) =z|j. (23) 

A first remark in this direction is that, if / = [n] is not a FSC for z, we can at 
least start the avalanche in the definition of recurrent by identity test, i.e. we have 
at least one site i — a{l) which is unstable. Indeed, we have at least one i such that 
z^ > /max(N)«- Then, by 1^, (z + b+), > /max([n]), + b+ = z^. 

As holds for any set /, the reasoning above works for any restricted sandpile 
A\i. Suppose to have a configuration z^^^ which is recurrent in A\i by absence of 
FSC's. Then we have at least one site i € I which is unstable, because we have at 
least one i such that > fmi,AI)i^ and {z^''> + b+{I))i > fraa,^{I)i + = Zi. 

This allows to construct an induction. Let /' = I \ i, v^'^ = (tiZ,'^-'^)\i' and 
z^'"> = n(y(^">). For j e /', vf'^ = z^p - A,j > z^ . Remark that the definition 

of fmax(>/) is the same on any restricted sandpile Ai with I ^ J. Thus, as z^^^ is 
recurrent for absence of FSC's, and this property is preserved under relaxation (by 
Proposition [T]) , also z^^ ^ is recurrent by absence of FSC's, on A\i'. This, together 
with equation (1^ . gives the induction step. 

In summary, we can perform the complete avalanche tcr(n) ■ • • icr(i)j algorithmically, 
by initialising z^^^ = z + b"*" and Iq = [n], and, for a = 0, . . . , n — 1, is any site i at 
which z^"^ is unstable (which is proven to exist by the reasoning above), la+i — /a\*o, 
and z('^+i) = 7^(^i^z('^)). 

This completes the equivalence of our three definitions of recurrent configurations, 
thus from now on we will omit to specify the defining property. 

A directed graph can be associated to a sandpile, such that ~Aij directed edges 
connect i to j. This completely encodes the off-diagonal part of A. The remaining 
parameters, in particular b"*" and b~, can be encoded through directed edges incoming 
from, our outgoing to, a special sink vertex. In the undirected case, A = A'^ , a 
bijection, called Burning Test, relates stable recurrent configurations to spanning trees 
[7]. A generalisation, called Script Algorithm extends this result also to directed 
graphs even in presence of greedy sites. Through Kirchhoff Matrix- Tree Theorem, one 
thus gets in this case that the number of stable recurrent configurations is det A. 

The configuration p = z-l-1 — 7?,(z-|-l) has two interesting properties: Pi > 1 for all 
i, and p ~ 0. This implies that, for every z S Z", the iteration of the map z TZ{z+p) 
must reach a fixed point, in R, the subset of S containing recurrent configurations 
(this would be true also with b+ instead of p, but slightly harder to prove). Indeed, 
calling c = maxi(zi — + 1), as TZ{TZ{- ■ ■ TZ{z + p) • ■ • + p) + p) = Ti-iz + cp), and 
z + cp is unstable at all sites, TZ{z + cp) is both stable and recurrent (by toppling 
covering), so it must be in R. 

This reasoning proves that each equivalence class has at least one representative 
in R. In the case b"*" S fi, as both the cardinality of R and the number of classes are 
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det A, each equivalence class must have a unique representative in R. In particular, 
the representative in i? of is called the recurrent identity, and we denote it with 
the symbol Id^. This configuration can be found as the fixed point of the map z — >■ 
TZ{z + b+), started from [37138) . or of the map z TZ{z + p), or, under the mild 
assumption that z G ^2, more directly, with no need of iterations, by the relation [33] 

Id^ = 7^(z+ (z-7^(z + z))). (24) 

Indeed, as z — TZ{z + z) G J?, z + (z — TZ(z + z)) is recurrent, thus its relaxation is 
in R. 

In the z — gauge, the operation u © v := TZ{u + v) sends S x S into S, and thus 
defines a semigroup on this space. Furthermore, this operation also sends S x R ^ R, 
R X S ^ R, and, as a corollary, R x R ^ R. 

The charges behave linearly under this operation: Q(u© v) — Q(u) — Q(v) e Z". 
The unicity of representatives in R of the equivalence classes allows then to construct 
inverses, for the action on this space, and thus to promote to a group action on 
i?, of which Idr is the group identity. This structure was first introduced by Creutz 
[37138) . and then investigated in several papers [20I35I39I40] . 

The covariance of the notations allows to define the operation © in general. We 
have 

u © V TZ{u + V - z) ; u ©^l^ v := Te'i' (u + v - z) . (25) 

Note however that now the charges behave in an affine way, and only the translated 
versions, Q'(u) = Q(u) — Q(z) and Q"(u) = Q(u) — Q(z) respectively for the two 
operations, have no offset. In particular, the two groups induced by ffi and ©^, as well 
as the two sets R and R\ are isomorphic but not element- wise coincident, as the only 
natural bijection among the two makes use of the involution t. 

Call A4 = M[a-,al] the transition monoid generated by the a^'s and al's acting on 
our set of configurations. A generic element in A4 has the form 

A = ai--- ai a.2 • • • a ^ at • ■ • at ■ • • a 2k ■ ■ ■ a 2^ ■ (26) 

'l '«(!) 'l '^(2) 'l V(3) 'l '«(2fc) 



3 Multitopplings in Abelian Sandpiles 

Consider a situation not dissimilar from the BTW sandpile, with height variables Zi at 
the sites of a square lattice, but let the toppling ti occur if the gradient passes a given 
threshold. For example, choose to perform a toppling at j if « is a local maximum. 

Despite the simplicity of such a rule, abelianity is lost, as is easily verified. For 
example, on a configuration z as below, the two sites i and j are unstable, but tiZ 
and tjZ differ, and it can well be that both have no unstable sites. 
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We describe now another variant of the BTW sandpile in which abelianity is lost. 
We come back to toppling criteria based on having heights above given thresholds. 
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as in the ordinary BTW (and not gradients, as above). Now the rule is that, if two 
neighbouring sites have both height 3 or higher, both their heights are decreased by 
3, and the height of each of the 6 neighbours is increased by 1. Again, it is not difficult 
to imagine situations in which two sites i and j are unstable, but t^z and tjZ differ, 
and both have no unstable sites, as in the configuration below. 
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In this section we will describe a family of sandpiles for which abelianity is preserved. 
In these models we have toppling rules based on "heights above given thresholds", 
just as in the ordinary ASM and in the example above. In fact, these models generalise 
the ordinary ASM, because we have simultaneously all the ordinary toppling moves, 
and a (possibly empty) collection of multitoppling moves satisfying certain criteria. 

As we will see later on, the moves described in the second example would be a valid 
choice - and remark indeed how, including also ordinary topplings, the configurations 
on the left and right of the figure above are still unstable on j and z, respectively, and 
become identical after these topplings. 

Consider an ordinary ASM ^(Z\,z, z). For simplicity of notations, we set in the 
positive-cone gauge, so that site i is unstable if Zi > An. 

For a non-empty set / C [n], call A/ — ^-^j A^. A multitoppling operator tj 
can be associated to a set /. First of all, z is unstable for toppling / if Zj > {Aj)j 
for all j (z I (note how, in the positive-cone gauge, this coincides with the ordinary 
definition when |/| = 1). Then, if the configuration is unstable, it is legitimate to 
perform the toppUng tjz = z A/. Also note that the absence of greedy sites implies 
that, for all j S /, {Ai)j > 0, while for all j ^ I, {Ai)j < 0. 

Consider a collection C of non-empty subsets of [n] . The interest in multitoppling 
rules for the Abelian Sandpile Model is in the following fact 

Proposition 2 Suppose that, for every I, J G C, the sets I' — I \ J and J' ~ J \ I 
are either empty or in C. Then, if z is unstable for both I and J, tjz is unstable for 
J' and tjz is unstable for I' . 

As clearly tjdjz = tjrtjz, it easily follows 

Corollary 2 In the conditions of Proposition\^ the operator TZ is unambiguous. 

Antitopplings t\ and antirelaxation TZ^ are defined just as in the ordinary case, e.g. 
through the involution i, which is still defined as in (|5]). A configuration is stable if 
no toppling or antitoppling can occur (this coincides with the definition of S in the 
ordinary case). Note that, if C does not contain the atomic set {i}, S is either empty 
or of infinite cardinality (because, if z e S*, also z + c Si G S for any c € Z). In order 
to exclude this pathological case, we will assume in the following that C includes the 
set Co = '^^ atomic subsets, i.e. single-site topplings. In this case, as 

stability w.r.t. £ is a more severe requirement than stability w.r.t. C C £, we have 
that is a subset of the set Sq of stable configurations in the associated sandpile 
with only single-site topplings, and thus of finite cardinality. As we will see later on, 
for any set C as in Proposition [5] the set S is non-empty, and actually contains a set 
isomorphic to R, thus it has cardinality bounded below by |i?o| — det A, an above by 

i^oi-n,(^.-^i + i)- 

If we require both that £ 3 £o, and satisfies the hypotheses of Proposition [31 we 
get that £ is a down set in the lattice of subsets, that is, for all / G £ and H C C 
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non-empty, also H E C (this is trivially seen: with notations as in Proposition [21 take 
J and J = / \ H). 

For a multitoppling sandpile A ~ A{C), we will call — -^{^o) the associated 
single-site toppling sandpile. 

The concept of recurrent configuration has to be reanalysed in this context. The 
various alternate definitions are modified (in the more complicated situation, but 
under the simplification of the choice of positive-cone gauge) into 

Definition 4 A configuration z is recurrent by identity test if there exists an ordered 
sequence (/i , . . . , /fc) of subsets of [n] , constituting a partition of [n] ( i. e. for all i G [n] 
there exists a unique a such that i € la), such that tj^ ■ ■ ■ i/^t/j (z -I- b+) — z, is a valid 
toppling sequence. 

Definition 5 A configuration z is recurrent by toppling covering if there exists a 
configuration u, such that z — tj^ ■ ■ - tj-^u is a valid toppling sequence, and each site 
i is contained in at least one of the la 's. 

Definition 6 A configuration z is recurrent by absence of FSC's if for every set I 
of sites, there exists J E C with L = I (1 J =/= , and Zj > X]ie/-^L "^^^ ^ ^ ^■ 

All the reasonings are the immediate generalisation of the ones already done in Section 
12.31 We just report here the appropriate modifications in Proposition [T] (recall that 
in our gauge SI = S'^, and is left stable by the topplings). 

Proposition 3 Let u G i? and v = tj^t]^ -^ ■ ■ ■ t/^u. Define A ~ Ui<a<*: ^a- For any 
set B, there exists a non-empty set H Q B, such that B C_ H U A, and, for all s E H, 
> ^J2ie{AnB)^H Dis- 



proof. We wiU produce explicitly a valid choice of H . U B ^ A, let H = B \ A 
and u' = u. Note that {B n A) \ H = B n A in this case, li B D A, for aU i E B, 
call t(z) the maximum 1 < a < A: such that i E la, then call r = maxiT(i), and 
J — lT(^i-^, i.e., the multitoppling that covered any portion of B more far in the past. 
Call u' = tj^ ■ ■ -t/jU, the configuration obtained after this last multitoppling. Note 
that all the entries of u' are non-negative. In the remaining part of the avalanche, 
for some non-empty set H C no more topplings occur. Conversely, all the sites 
j E (AnB) \ H do topple at least once (possibly in a multitoppling event). Thus we 
have, for each s E H, 



> - ^ A. 

ie(snA)\// 



(27) 



as was to be proven. □ 
We should modify the concept of forbidden sub-configuration along the same lines. 

We want to produce pairs (/, f) such that, if z is recurrent by toppling covering, 

z\i ^ f. For a given I, a vector f has the property above if, for all J G £, J C I, 

f|,7 ^ ^J2k£i-^j^k\,j- Note that the fact that £ is a down set has been used to 

restrict the set of J's to analyse. 

A collection of the forbidden pairs (/, f ) with smallest |/|, such that no other f 

exists with i' y i and (/, f ) a forbidden pair, in the BTW sandpile with multitoppling 

rules on all pairs of adjacent sites, is as follows: 
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Note that, at difference with the single-site sandpile, in general there is not a unique 
fmax(^) such that (/,f) is a forbidden pair if and only if f ^ fmax(^)- 

At the level of the monoid 7W[aj, aj], we have some extra relations associated to 
multitoppling rules. For example, consider the action on S'^. While in the single-site 
sandpile we have relations ([Ti)) . in the multitoppling case we also have relations of 
the form, for each / e 

n«r^'=n«7'^'^^- (28) 

Note that multiplying left-hand and right-hand sides of relations (ITi)) for all j € I we 
would have got a weaker relation, in which the two sides of (j28p are multiplied by 

n «7^*^ ' (29) 

(recall that at the level of the monoid it is not legitimate to take inverses of a^'s, and 
that, in the equations above, all the exponents are indeed non- negative). 

At the level of the equivalence relation ~, and thus of charges Q(z), nothing 
changes. In particular, it is easy to see that u ~ v in ^ if and only if u ~ v in ^o- 
This also leads to the fact that there is exactly one stable recurrent configuration 
per equivalence class, and that © defines a group structure over R, just as in the 
single-site toppling sandpile Aq = A{Co) associated to the multitoppling sandpile 
A^A{C). 

Note however that the set R is different from the set Rq of stable recurrent config- 
urations for Aq, and the fact that they have the same cardinality results from subtle 
compensations between stable/unstable and recurrent/transient configurations. For 
example, in the BTW sandpile with multitoppling rules on all pairs of adjacent sites, 
the set R loses the configurations with adjacent pairs | 3 | 3 | , which are now unstable, 
but gains the configurations with adjacent pairs | o | o L which are not FSC's anymore. 

A natural bijection between R and Rq, preserving the group structure, is obtained 
by associating to z G i?o the configuration TZ{z) G R, where TZ is the complete 
(multitoppling) relaxation. 

In the case of a tight sandpile, the set containing the conjugate of the stable 
recurrent configurations in the single-toppling sandpile coincides with the set of 
recurrent configurations of the sandpile in which C contains all the subsets of [n] . In 
this case, there are no stable transient configurations, and the condition for z being 
stable w.r.t. any toppling /, i.e. z|/ ^ A/ for all non-empty / C [n], is related, 
by conjugation, to the condition of not having FSC's in i-e., (tz)|/ fmax(-^) 
(because fmax(-f) = (z — A/)|/, and, in applying the definition of t, we should recall 
that the multitoppling sandpile is formulated in the positive-cone gauge). 

Now consider an anomalous relaxation process, p/, which may perform a multitop- 
pling rule / only at the first step, if possible, and then perform a single-toppling 
relaxation with TZq. Such a process is unambiguous, but different processes may not 
commute, i.e. pi{pj{z)) ^ pj{pj{z)) in general. 

Nonetheless, take a whatever semi-infinite sequence (/i, I2, Is, ■ ■ ■) of elements in 
C \ Co, such that all elements in C \ Cq, occur infinitely- many times. It is easy to 
see that, for all z, there exists a truncation time t — t{z) such that, for all s > t, 
Pi^Pi^-i ■ ■ ■ Pii (z) = 7?.(z), and in particular piTZ{z) = TZ{z) for all / G £. 

The interest in these processes pi is in the fact that their action is strongly related 
to the action of the idempotents a-a^. Before stating and proving this relation in 
precise terms, it is instructive to investigate first how this works in the case of the 
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BTW model. One easily recognizes that aja^z = z if < 3, or = 3 and Zj < 

3 for all the neighbours j of i. In the first case, no topplings or antitopplings are 
involved, while in the second case exactly one toppling and one antitoppling at i occur. 
Conversely, if Zi = 3 and Zj = 3 for some neighbour j, causes an avalanche for 
which a valid sequence of topplings may start with (i, j, . . .), i.e. a^z = tj^ ■ ■ • U^tjtiZ 
for some (13, . . . , z^). The effect of the two initial topplings is identical to the effect 
of a niultitoppling at a pair And, crucially, also the if condition coincides with 
the one for a configuration to be unstable w.r.t. the multitoppling at {i,j}. Thus, 
a configuration z G is left stable by the application of a|a,- if and only if it is 
stable w.r.t. both the toppling i and the multitopplings {i,j} for all j neighbours 
of i. This proves Theorem [5] in the case of the BTW sandpile, and characterizes 
y(z,/) as the result of the relaxation of z, in the multitoppling sandpile for which 
C\Co = {{i,j}}^^j^ d(^,,) = i ■ 

The proof in the general setting, that we present below, is completely analogous. 
Proof of Theorem [2J One finds that ala^z = z if < z^, or Zi = and zj < 
'Zj+Aij for all j ^ i- Again, in the first case no topplings or antitopplings are involved, 
while in the second case exactly one toppling and one antitoppling at i occur. 

Conversely, if the conditions above are violated, Oi causes an avalanche for which 
a valid sequence of topplings may start with (i, j, . . .), and the effect of the two initial 
topplings is identical to the effect of a multitoppling at a pair 

The if condition for the multitoppling {«, j} to occur is that Zi > 2.^ + Aji and 
Zj > 'Zj + Aij. The condition for the avalanche to involve topplings at i and j is 
Zi = 'Zi and Zj > 'Zj + Aij . These two sets of conditions are certainly simultaneously 
not satisfied if Aij — 0, thus we can restrict our attention to sites j such that Aij < 0. 
In this case, the two sets coincide if and only if Aji = — 1. Thus, in order to make 
the sets coincide for all sites i, we need that A is symmetric, and all the non-zero 
off-diagonal entries are —1, as required in the theorem hypotheses. 

The configuration y(z, /) is the result of the relaxation of z, in the multitoppling 
sandpile for which C \ Cq ~ a --1 ' ^ 

4 Discussion of Markov Chain dynamics involving both sand 
addition and removal 

In this section we discuss several dynamics involving the operators and a]. Each 
example has a different theoretical and phenomenological motivation, and is intended 
to describe a different feature of out-of-equilibrium steady states. To keep the visual- 
isation simple, all our examples are variations of the BTW model, on portions of the 
square lattice and with heights in the range {0, 1, 2, 3}. 

4.1 A reminder of the Karmakar and Manna protocol 

The first and most natural dynamics in this family is the modification of the BTW 
sandpile, where, with probability p and 1 — p an or an a| move is performed, at a 
randomly chosen site i. This dynamics has been investigated in detail by Karmakar 
and Manna in ISTl. For p > 1/2 and p < 1/2, the system has features resembling the 
ones of the ordinary BTW model, and of its symmetric image under the involution 
Zi o 3 — Zi. For example, for p > 1/2 the avalanches have an algebraic-tail distribution 
(the exponent varies with p) , while the anti-avalanches have a distribution with most 
of the support on values of order 1. Of course, these two features are swapped for 
p < 1/2. In a window near p = pc — 1/2, for which several scaling exponents have 
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been investigated numerically, aspects of a new criticality phenomenon emerge. For 
example, at p — Pc the distribution of sizes of both avalanches and anti-avalanches 
seem to follow a stretched-exponential. 

As another example of new emerging features, both the variance and the correla- 
tion times of the total mass in the system increase when pc is approached, and the 
variation in p of the average total mass diverges at Pc, with a new scaling exponent. 

4.2 Dynamics preserving the total mass 

Our first dynamics is a minor modification of the one described above, in the micro- 
canonical ensemble with fixed total mass, at its critical value. Among other things, 
this allows to use periodic boundary conditions. Indeed, the dynamics studied by 
Karmakar and Manna, even at p = pc, ii implemented on a graph with no bound- 
ary, may produce an infinite avalanche or anti- avalanche, while in a dynamics with 
conserved mass, with randomized initialisation, if the initial density is in the range 
3 — < p < p*, where p* is the critical density (having value p* — 2.125288... 
[41142) . note that, for a long time, it was mistakenly believed that p* = 17/8), in the 
large volume limit we are protected from infinite avalanches. 

We want to emphasise here the behaviour of thermalisation process at short times, 
that shows coarsening, and an aggregation phenomenon, which is purely dynamical. 

More generally, under some aspects, our dynamics compares to the one of Kar- 
makar and Manna roughly as Ising at fixed magnetisation under Kawasaki dynamics 
relates to Ising heat-bath evolution. Recall however that the presence of conserved 
quantities has more dramatic effects in non-equilibrium systems than it has in Boltz- 
mann equilibrium theory (see e.g. the very different features of the Driven Lattice 
Gas, at conserved number of particles or with open boundaries [44j ) . 

The precise time evolution is as follows. At each time, one of the 2n^ operators 
{ajflj , a^Oj} is chosen uniformly at random, and applied to the previous state. For the 
examples shown in Fig. [TJ we have chosen as initialisation the one best approximating 
a uniform distribution of the mass, = 1 if a; -I- y is odd and 2 if it is even. As 

seen from the figure, the coarsening is already evident at the level of the height 
variables. It is however much more evident at the level of a local parameter suggested 
by the Burning Test. The ordinary Burning Test cannot be performed here, as we 
have no boundary. We choose as effective boundary the set of sites with maximal 
height. More precisely, we record the sites that perform at least one toppling, in the 
avalanche following the replacement 3 — > 4, and colour them in blue. Similarly, we 
record the sites that perform at least one anti-toppling, in the avalanche following 
the replacement ^> — 1, and colour them in red. In case both a toppling and an 
antitoppling have occurred, we break the tie by observing the initial height. In case 
neither a toppling nor an antitoppling have occurred, we use a neutral light-yellow 
colour. Thus, pictorially, blue, red and yellow regions correspond to regions which 
are recurrent, anti-recurrent, and simultaneously transient and anti-transient. The 
resulting analysis of the configurations is shown in Fig. [2] 

4.3 Dynamics with idempotent operators 

Our second dynamics starts from the maximally filled configuration, Zi — i for all i, 
and acts with the idempotent combinations a\a^ at randomly-chosen sites. 

As we know from Theorem[21 this dynamics is absorbent on a unique configuration, 
identified with the multitoppling relaxation of the initial configuration, where pairs 
of adjacent sites both with height 3 are unstable. 
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Fig. 1. Time evolution of the Markov Chain described in Section 14.21 Here L 



256, and times shown are t 
(red, orange, cyan, blue). 



L\AL\QL^ 



,U4:L\ Color code: (0,1,2,3) 




Fig. 2. The local parameter associated to burning-test defined in the text of configurations 
in Fig. [U 
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Fig. 3. Configurations obtained with the procedure described in Section [4.31 Top, middle 
and bottom row correspond to t = 32, t — 1024, and to the fixed point of the dynamics. The 
columns show different domains. From left to right: a circle; an ellipse with axes rotated by 
arctan 1/3 w.r.t. the cartesian axes, and ratio 2 between height and width; a typical algebraic 

7; a square. 



curve of degree 3, more precisely 4x^ + 4j/^ + Sxy 



■ 2x^y + 



Again, interesting features emerge at short times, when the configuration takes 
the form of a web of strings, satisfying a classification theorem and a collection of 
incidence rules P7] . 

On the other extremum of the dynamics, at the fixed point we have configura- 
tions showing remarkable regularities, in the form of patches, that is, a local two- 
dimensional periodicity on portions of the domain ^^51. When more patches are 
present, they follow an incidence rule [26]. 

If the initial domain is an elliptic portion of the square lattice, a specially higher 
regularity emerges. Say that the linear dimension of the domain is of order L, and 
the slope of the symmetry axes is a small rational p/q (with both p and q of 0(1) in 
L). Then, in the limit L ^ oo, we observe the emergence of a very simple structure 
of patches and strings: we have a unique patch, crossed by strings of a unique type, 
parallel to one of the two symmetry axes. This fact is in agreement with the general 
theory developed in |25I26I27] , as the toppling vector at a coarsened level is a quadratic 
form in the coordinates x and y, that should vanish at the boundary of the domain, 
and the contour lines of quadratic forms are conies, i.e. plane algebraic curves of 
degree 2. 

In Fig. [3] we present configurations obtained with the procedure described above, 
starting with the maximally- filled configuration, = 3 for all i, on portions of the 
square lattice of various shapes, in order to highlight the features outlined above: a 
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disk, an ellipse, a smooth domain which is not a conic (it is an algebraic curve of 
degree 3), and a square. 

4.4 A simple deterministic dynamics 




Fig. 4. Examples of configurations, for L = 60, observed in the dynamics described in 
Section|331 The domain is repeated periodically 2x2 times, in order to improve visualisation. 
The different rows correspond to densities p = f.5, f.75, 2, 2.fff.... The three left-most 
columns show, respectively, configurations at time t — 3000, 6000 and 9000, obtained by 
starting from periodic configurations prepared by using the tiles (12)) (21)' (2)' ^^'^ 
/ 3 2 2 N 

(222). The last column shows, for comparison, the configurations at t = 3000 obtained 
by starting from random configurations at the given densities. 

Our third and last example of dynamics is deterministic, and corresponds to the 
sand flow from a unique source site i to a unique sink site j in the lattice, through 
the iterated application of the operator a^^a^ . 

Besides producing beautiful drawings, this dynamics has an interesting phenomeno- 
logical feature. Contrarily to the nai've intuition of a smooth hydrodynamic flow from 
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the source to the sink, in typical configurations at the steady state, the density pro- 
file presents a shock, between a spatial region analogous to ordinary BTW and a 
region analogous to the involution image of the BTW, both regions presenting large 
avalanches (of the appropriate type) . 

This phenomenology seems an analogue of the shock profiles observed in various 
non-equilibrium systems, most notably a two-dimensional analogue of a feature of the 
TASEP 45i, an exactly-solvable model in one dimension, in the case in which two 
species of particles are considered [15] . 

We note that, in comparison with the previous realisations, this dynamics presents 
short thermalisation times and very moderate aging, (both the thermalisation time 
and the aging effect seem to become more relevant as p approaches p*). A partial 
justification is the fact that this dynamics is deterministic and thus, after a transient, 
it follows a periodic orbit, although the growth of the period with the system size is 
extremely fast. Recall that the length of the orbit of the operator aj^at, acting in the 
torus G = X X • • • X Z^^ of recurrent configurations, for the abelian version of 

the model, is a divisor of the length of the orbit of a]aj acting on S, due to the fact 
that the equivalence classes are preserved in the framework with both topplings and 
antitopplings. 

We presents configurations realised on a L x L torus, with sink and source at 
antipodal points, i.e. acting with the operator ajoj = 0(^/2 l/2)'^(o,o)i which is the 
simplest and most symmetric realisation of the ideas above. We studied both the case 
in which the initial state is periodic, and the case of random initialisation with given 
prescribed density. The configurations are presented in Fig. U) 
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